Distribution of the resistance of nanowires with strong impurities 
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Motivated by recent experiments on nanowires and carbon nanotubes, we study theoretically 
the effect of strong, point-like impurities on the linear electrical resistance R of finite length quan- 
tum wires. Charge transport is limited by Coulomb blockade and cotunneling. Ini? is slowly 
self-averaging and non Gaussian. Its distribution is Gumbel with finite-size corrections which we 
compute. At low temperature, the distribution is similar to the variable range hopping (VRH) 
behaviour found long ago in doped semiconductors. We show that a result by Raikh and Ruzin does 
not apply. The finite-size corrections decay with the length L like 1/lnL. At higher temperatures, 
this regime is replaced by new laws and the shape of the finite-size corrections changes strongly: if 
the electrons interact weakly, the corrections vanish already for wires with a few tens impurities. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.20.Ee 
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Understanding the charge transport in materials with 
impurities has been a long standing problem which at- 
tracts new interest because of recent experiments on sys- 
tems of nanoscopic size III] . In the past, the case of doped 
semiconductors has been extensively studied [5] . Because 
of impurities, the electrons are, in equilibrium, in local- 
ized states. These states are randomly spread in energy 
and space. At low temperature, conduction results of a 
sequence of thermally activated hops of the charge car- 
riers, which tunnel from one state to another, borrowing 
the energy difference from the phonons bath if needed. 
The spatial range of these hops results from a balance 
between tunneling (favoring short ranges) and activation 
(finding a state close to the Fermi level is more likely if 
the range is long). As a consequence, the low-field re- 
sistance R of amorphous semiconductors follow Mott's 
VRH law, Ini? ~ in dimension D, with pe- 

culiarities in the special case of ID 0, In presence 
of long-range Coulomb repulsion, as for crystalline semi- 
conductors, the states close to the Fermi surface are rare 
and this law must be replaced at very low temperatures 
by the Shklovskii-Efros VRH law, Ini? - T^^/^ 

Here we deal with ID systems of nanoscopic cross- 
sections (nanowires) with point-like, strong impurities, 
and repulsive interactions between the electrons, at low 
applied electric field. In a pure nanowire, or between two 
impurities, the physics is described by the Tomonaga- 
Luttinger liquids (LL) model rather than a Fermi liquid 
because, in a ID fermionic system, the excitations have 
a bosonic character. The case of a few strong impuri- 
ties, or of a uniform density of weak impurities, have 
been studied 0, Q. R was found to behave like pow- 
ers of T at least in some regimes, with exponents which 
depend on the interaction parameter K of the LL. As 
was recently shown Q , in the case of many strong impu- 
rities. Coulomb blockade [l[ and statistical effects yield 
different laws: at low temperature, a VRH-like regime 
is predicted, while at medium temperature a new, more 
complicated law (which can, on a small range of T, look 



like a power law) should replace it. 

In this Letter, we first confirm numerically the exis- 
tence of these two regimes for the average of In i?. But 
they might be hard to distinguish due to huge sample to 
sample fluctuations. Therefore, we study the full distri- 
bution of Ini?. R and Ini? are very slowly self-averaging 
and non Gaussian. Unless the wires are very long, the 
distribution of In R is Gumbel with corrections that decay 
like 1/lnL at low T. They are much weaker at medium 
T and have a different shape. We think this will help 
uncover experimentally the two regimes. Finally, we dis- 
cuss the ergodicity hypothesis often used in experimental 
studies. 

Our plan is the following: after introducing the model 
and the simulation technique, we show numerically that 
each sample has a high T regime where its resistance is 
proportional to the length and a "non extensive" regime 
where the voltage drop occurs essentially on a few breaks. 
We check the validity of a percolation-like approximation 
to R. Then we give our conclusions for the average of In R 
and derive the distribution of In R. 

Model. Our model is similar to the one of Refs. [1,0]. 
We consider spinless electrons. Their interactions are 
short range because of screening by a nearby gate, which 
gives a finite capacitance per unit length C to the 
nanowires. Each wire is cut into a series of weakly cou- 
pled quantum dots by TV impurities. They are placed 
randomly according to a Poisson process with mean spac- 
ing I large compared to the Fermi wavelength \p. In the 
quasiclassical limit, which is relevant because the impu- 
rities are strong, each dot j can only have an integer 
number, g^, of electrons. Thus the charge transport pro- 
ceeds by discrete hops, like in disordered semiconductors, 
and is limited by activation and tunneling. But unlike 
the localized states in the VRH model, here the ener- 
gies of the states are correlated with their positions. In- 
deed, in the ground state, qj is the closest integer q-p 
to the mean charge Qj = Ij/Xp + ji/Aj, where fj, is 
the chemical potential and Aj = l/{Clj) is the charg- 
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ing energy (we use e = fee = !)■ The energy of dot j is 
^jilj) = i'lj~Qj)'^^j/'^ 0- The typical charging energy 
A := 1/(IC) is the natural energy scale of the problem. 
In addition, each charged excitation qj ^ gj"'' may itself 
be neutrally excitated by an energy difference ~ KAj. 

Having in mind experiments where a four-point mea- 
surement is performed [ij, we do not pay attention to 
the contacts between the wire and the two leads. We 
assume for convenience that each impurity has the bare 
transparency l/(e7r) « 0.12. According to [oL [ill ll^. 
because of the LL effects, its effective tunneling trans- 
parency may be found by an instanton calculation, yield- 
ing exp(-Sj,j+i) with Sjj+i = 1 + {K^^ - l)(crj-i + aj) 
with aj = ln[K/{CXF)/ max{KAj,T)]. The tunneling 
process is limited by neutral excitations at low T and 
by thermal fluctuations at high T. This extends readily 
to the probability exp(— Sj^fc) of the cotunneling of k — j 
electrons, each through one of the impurities j, . . . , k — 1 
(which is formally equivalent to the hop of one electron 
from the dot j to the dot k). At low temperatures, co- 
tunneling allows to shortcut the narrow dots with high 
charging energy — this is why In R increases slower than 
l/T as T goes to 0. 

Putting things together, we find the probability that 
an electron hops from dot j at Xj with initial charge qj 
to dot k at Xk with initial charge qk under field F: 



''-foE[exp{E/kBT)-l]- 



where E = E^iq^ - 1) - Ej{qj) + Ek{qk + 1) - Ek{qk) - 
eF[xk — Xj). The energy difference ^ is pro- 
vided/absorbed by the phonons bath; the field F favors 
one hopping direction. 70 depends on the coupling with 
phonons [l3|. 

In a steady regime, the average occupation probabil- 
ity of each state should be constant, therefore we write 
down for each state a balance equation between incom- 
ing and outgoing charge flows. Solving this in a general 
out-of-equilibrium setting is difficult, but, following the 
idea of Miller and Abrahams [l3| , we restrict to the small 
F regime and simplify these equations to the first order 
in F. The resulting equations are the Kirchhoff's laws 
for a resistor network between nodes labeled by j, qj . In 
the following we are concerned with the low temperature 
regime T < lOA and we assume that only the first two 
charged excitations on each dot play a role. If this may 
not be true for all dots, it will certainly be the case for the 
dots which contribute the most to the resistance by their 
high charging energy. Moreover, we further simplify p?] 
to a network with one node per dot j and resistances Rj^k 
equal to the minimum over qj and qk of the resistances 
between the states j, qj and k^q^. In this network, each 
node is connected to all others, although connections be- 
tween remote dots are likely to be very resistive. 

Simulations. We performed simulations by drawing at 
random the positions of the impurities of 10000 wires for 



each value of impurities number TV and mean spac- 
ing I. The low field resistance R of each wire is com- 
puted as the resistance of the equivalent network defined 
above; more precisely, we took care to manipulate Ini? 
rather than R to prevent overflows at low T where R 
become huge; arbitrary precision arithmetics was used 
when needed to extract subdominant values. Simple lin- 
ear solving of the Kirchhoff's equations fails at low tem- 
peratures because of numerical instabilities, but repeat- 
edly eliminating one node of the network by applying 
a generalized star-triangle transformation 1^ until only 
two nodes are left is efficient and stable. 

For each wire, we also computed the resistance of the 
"best path", R^p, which is the less resistive of the sub- 
networks without loops connecting the nodes and 1 
and going always in the same direction, in analogy with 
the percolation approach to VRH . A link in the best 
path is termed a hop. We found that In i?bp is a very 
good approximation of Ini? up to T ~ lOA. Clearly, 
Rhp > R because the full network has derivations around 
the best path. Empirically, In i?f,p/ In i? < 1.01 always, 
and this ratio goes to 1 as T decreases. On our curves 
the averages of hiRbp and Ini? can't be distinguished. 
In the following we use Rbp for our analytical discussion 
because it is a sum of individual resistances. 

From the simulations we conclude first that each sam- 
ple has two well separated regimes: at high T the voltage 
(1) drop is almost uniform and the resistance Rbp « i? is ex- 
tensive (proportional to the length), but at low T the 
voltage drops mainly in one or a few regions, the breaks 
(Fig. [1]) , which makes R non extensive unless the wire is 
extremely long [l3| . The positions of the breaks strongly 
depend on T and their width increases as T goes to zero 
and cotunneling through more and more impurities takes 
place. Fig. [1] shows this for one particular sample but it 
is true for all but pathological samples and for all values 
of K. The very low temperature regime T < O.OIA is 
cut into sample-specific ranges where R{T) is activated, 
as was observed experimentally [ij : in each range of tem- 
perature, a different break fixes the voltage drop. 

The second conclusion is a confirmation of the results 
of Q for the average over many samples (In R) : a VRH- 
like law (Ini?) ^ 1/Vt holds at low T but has to be re- 
placed by a new regime at intermediate T where l/(ln_R) 
is a afHne in InT (Fig. [2]). However, unless the wire is 
very long, the temperature range where this new, non 
VRH-like law can be observed is small and it is in the 
cross-over between low and high temperatures (cotun- 
neling and sequential tunneling), so that distinguishing 
the two laws is hard, all the more than the sample to 
sample fluctuations are large (Fig. [2]). But we shall see 
that looking at the full distribution of Ini? may help. 

Analytical results for the distribution of In R. From 
now on, we restrict to the low T regime, where breaks are 
observed. Taking for simplicity uniform s and Xp = 0, 
and neglecting correlations between hops, the probability 
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FIG. 1: The gradient of the electrical potential V{x) along 
the best path in one random sample with A'' = 100 impurities, 
interaction parameter A" = 0.3 and mean impurities spacing 
7 = IO'^Af. For each value of the temperature T, we apply a 
voltage unity to the wire and we compute the gradient V' of 
V on each hop along the best path (see text). Then we color 
this hop with a grey intensity depending on V' (the bar above 
the figure indicates the intensity coding; white is null gradient 
and black is the maximal gradient over the whole figure) . Top 
inset: The sum of the squares of the voltage drop of each 
hop along the best path is 1 at low temperature (one hops 
dominates) and at high temperature (A'^ hops with drop 
1/N). For each of 10* samples from the same distribution, 
we look for the temperature T* where it is equal to l/\/7<f. 
This separates the non-extensive and extensive R regimes (see 
text). Strong lines: median value of T* for several A'" values 
(solid: K = 0.1, dashed: K = 0.7). Dotted lines: boundaries 
of the 68% confidence interval on T* at a given A'^. Bottom 
inset: Plot of In R w.r.t. 1/T. The curve is made of activated 
segments. In each one, the voltage drop occurs in one break 
which is identified by a letter on the main plot. 



that a given hop on the best path has In i? < u is 



P{i?< e''} = 1 



k 



T 



(In u— J s)] with fc = [_u/s\ (2) 



where f{x) := 1 — x [1 — exp(— 1/a;)]. The number 
n of hops in a sample is proportional to N (at low T, 
n w y^T/AiV) and, because of the equivalence of sta- 
tistical ensembles, we can replace n or by i in our 
asymptotic estimations. The most resistive hop has typ- 
ically a resistance 7?max such that P{R < i?max}" ~ 1/2. 
As long as i?max is large compared to the average hop 
resistance times n we are in the non extensive R 
regime, Rbp ~ i?max — this regime extends at least to 
wires with several tens of thousands of impurities. If 
Rbp/s > 6, — In P[i? > e"] is well approximated by an 
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FIG. 2: Solid line: numerical results for the inverse of the 
average of Ini? over 10* wires similar to the one of Fig. [1] 
as a function of T/A. The error bars are smaller than the 
line's thickness. The thin, dotted lines above and below in- 
dicate the boundaries of the 68% confidence interval at each 
temperature. Strong short-dashed line: high T asymptote 
and upper bond, reached when the electrons tunnel sequen- 
tially through all impurities. Strong long-dashed line: low T 
asymptote where In(ii) follows the VRH law, as found from 
the inset. The thick straight line indicates the validity of the 
new law predicted in [3| in some range of T . Inset. Dots: 
plot of the average In R w.r.t. ^J'KJt to show that, at low T, 
the VRH law is obeyed. The error bars are smaller than the 
symbols' size. Straight line: linear fit , which is used on the 
main plot as VRH asymptote. 



integral, which itself takes two forms Q: (a) Tu^/(sA) if 
mT/A < 1 and (b) (m/s) ln(uT/A) if uT/A > 1. The 
conditions Rbp/s ^ 1 and T/A In i?bp <C 1 are always 
fulfilled at (very) low temperatures, say T < O.OOlA. 
Furthermore, we verified numerically that the cumula- 
tive distribution function (CDF) of In R of the hops along 
the best path of the wires (with nonuniform s) is in very 
good agreement with the shape exp(— u^) once rescaled. 
Then In i? w In Rbp has for large n a Gumbel distribu- 
tion because it is the maximum of n random variables 
with exponential decay of the probability density 17, 1^. 
Its average behaves like y/s{lnn)A/T and its standard 
deviation a like a/2/ (In n)T/{As) where a := 7r/V6, 
which is similar to the results of Lee et al. for the VRH 
model [4j and agrees with the low T result of 0] . Both R 
and In R are non gaussian but self-averaging. We prefer 
to study Ini? because it is easier to sample numerically. 
We computed the CDF of the centered and reduced In R 
with finite n corrections 11911: 



Ini? - (Ini?) / M n 

P{ ^ <x}^ exp[-5r(a;)]{l - 



(3) 



g(x)[7rV24(l - x^) + a;C(3)/(2p)]/ Inn + 0[l/{lnnf]} 
where g{x) = exp(— 7 — ax), C(3) ~ 1.202 and 7 ~ 0.577. 
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FIG. 3: Solid line: CDF_of InR for wires with TV = 100 
impurities, K = 0.7 and l/Xp = 1000 at temperature T — 
0.1 A from our simulation. Dotted line: the same for one 
sample but /j, varying uniformly in the rage — 1..1 to show 
that the ergodic hypothesis is only approximate. Long/short 
dashed line: Gumbel/Gauss distributions having the same 
average and variance than this CDF. The Gumbel distribution 
is almost superimposed on the data. Insets. Solid line: finite- 
size effects for the same wires but K — 0.3 and at T = 0.1 A 
resp. O.OOIA (difference between the CDF of the centered 
and reduced Ini?, and a Gumbel distribution). Long dashes: 
eq. (|3} with n — 100. Dots: shape of the finite-size effects 
after Raikh and Ruzin's result. 



This expression is successfully compared to the numerical 
CDF of In R in the inset of Fig. [31 Since 1/ In n goes 
slowly to zero, the actual Ini? of each sample will likely 
lie away from the average (In R) even if n gets large. This 
is why activated segments are still visible on long wires 



Pursuing the work of Lee et ai, Raikh and Ruzin 15 1 
evaluated the distribution of Ini? in the case of ID VRH. 
Interestingly, their result tends to a Gumbel law for very 
short wires or at very low T. But the finite-size effects 
are in and do not match eq. nor our numerical 

results (bottom inset of Fig. [3]). 

At not so low temperatures, say T > O.OlA, the 
CDF of individual hops doesn't agree with form (a) 
above but, if interactions are not too strong, with form 
(b). The VRH-hke law (Ini?) ~ l/VT is replaced by 
(In R) w s\nn/ ln[(sT/ A) In n] where the number of hops 
n is proportional to L as before. The difference between 
the averages in the two regimes is small compared to the 
fluctuations a ^ 1/lnlnZ/. But a strong change in the 
finite-size effects of the CDF of In i?, which is Gumbel- 
like, should be easy to detect: their amplitude is reduced 
by ln[(sT/A) Inn] i.e. more than 3 in our data and their 
shape is changed. Fig. [3] shows that, for K = 0.7 and 
only N = 100 impurities one can't tell the difference be- 
tween the CDF of In R and a Gumbel CDF. For stronger 
interactions (smaller K), say K = 0.3, a typical break 
has only a few impurities, (b) is not valid, and the CDF 



is hard to characterize, but the deviations from a Gumbel 
law can't be confused with the formula ([3]) of the VRH 
case (top inset of Fig. [3]) and might still be used as an 
evidence of strong impurities in a LL wire. 

In an experimental study of the distribution of In R, it 
is common to measure R for the same sample but several 
values of the gate potential: this changes fi, and the ac- 
tivation energies strongly depend on /i. This was argued 
to yield the same results as using many samples (see 
and references therein). However, in addition to vary- 
ing fi, we recommend to use as many samples as possible 
because we found numerically that this ergodicity hy- 
pothesis is only approximate and may lead to systematic 
errors (Fig. [3]). 
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